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Abstract 

The nonlinear evolution equation for the scattering amplitude of colour dipole off the 
heavy nucleus is solved in the double logarithmic approximation. It is found that if the 
Qh| initial parton density in a nucleus is smaller then some critical value, then the scattering 

amplitude is a function of one scaling variable inside the saturation region, whereas if it 
is greater then the critical value, then the scaling behaviour breaks down. Dependence of 
the saturation scale on the number of nucleons is discussed as well. 
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1 Introduction 



In the two previous publications |l], [| we solved the nonlinear evolution equation for high parton 
density QCD[[| |, |5|, ||, [7|, ||, [9|, [10], [U], [L2|, [L3| in the double logarithmic approximation and 



discussed its properties. We showed that at high energies there is a momentum scale Ql(y, A), 
where y = \b.(1/xb) is a rapidity, at which the parton evolution can no longer be described by 
linear evolution equation. This so-called saturation scale, has the following dependence on the 
energy and the number of nucleons (for DIS on a nucleus): Q 2 s (y,A) ~ Ase 4oisV , which holds 
in the double logarithmic approximation. In Ref. [p], |2|] we derived this result assuming that 
initial parton density is not very large. 

The A dependence of the saturation scale can be also obtained using simple qualitative 
arguments. Indeed, as one increases energy and/or decreases virtuality of the photon, the 
density of partons inside the nucleus grows until they begin to overlap in the transverse plane. 
The scale at which this occurs is proportional to Q 2 (y,A) and can be estimated from the 
equation 

Atix 2 xbG(xb, Q 2 ^ 



ttR a 



(1) 



where x is a colour dipole, xbG(xb, Q 2 s ) is the proton structure function and Ra is the nucleus 
radius. Since in the kinematical region intermediate between the linear and the saturation ones, 
the anomalous dimension of the gluon structure function 7 ~ ^|], TDj, [T3| and Ra ~ A 3, one 
obtains that Q 2 ~ A§. 

It was shown in Ref. |7| that the initial condition for the nonlinear evolution equation is 



given by the Glauber-Mueller formula] 15 



N(x,y = y ) = l-e-^ XBG ^^ s ^ , (2) 



3A 



where S(b t ) is a nucleus profile function and b t is an impact parameter. If b t < Ra, S 

— A 

(Woods-Saxon parameterization |]nj), which implies that initially the saturation scale Q 2 S is 



proportional to As . So, the question arises how it becomes proportional to As at large energy. 
There is an opinion that the only dimensionful scale in the high energy DIS is the scale set 
by initial condition. It would imply that the A dependence of the saturation scale set by 
initial condition Eq. (0) remains the same during the evolution since the evolution equation is 
conformally invariant. This however is not the case since there is another dimensionful scale 
which comes through the Glauber-Mueller initial condition Eq. (|2]). Indeed, in our estimations 
we neglected for simplicity the logarithmic dependence of the gluon structure function xbG of 
the nucleon on x 2 . However, this dependence on ln(x 2 A 2 ) introduces the second dimensionful 
parameter A which is the soft QCD scale. In other words, the relevant scale for the evolution 
in the nucleon is A 2 while that in nucleus is AjR\. The later one as we see is set by the nucleus 
size. The former one means that evolution can occur in a nucleon as well as in nucleus. 

The goal of our paper is to derive the saturation scale in the semiclassical double logarithmic 
approximation and discuss its dependence on A and energy. We consider the nonlinear evolu- 
tion equation derived by Yu. KovchegovQ, [7j in the framework of Mueller's dipole model [O 
and solve it using semiclassical methods (on characteristics) in the double logarithmic approxi- 
mation. We show that there exist a certain characteristic, which divides the whole kinematical 
region into two patches: one in which the evolution is linear and another in which the par- 
ton density saturates. The remarkable property of this characteristic is that the dipole number 
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Figure 1: The kinematical regions relevant in high 
parton density QCD. The initial condition Eq. (Q) 
for the nonlinear evolution equation Eq. (p|) is set on 
the line y = yo (dashed line). Solid curves 1,2 are 
the critical lines r = (\nQ 2 s /4A 2 ) for different nuclei: 
A\ < A 2 . Initial conditions for the solution to the 
left of the critical line are set along line 1 for nucleus 
A\ or along line ABC for nucleus A 2 . 



remains fixed along it. This characteristic coincides with the critical line in the double logarith- 
mic approximation^, [Tj|. Recall H that the critical line is defined by r = {ha.Q 2 {y, A)/4A 2 ), 
where r = (lnQ 2 /4A 2 ). 

In this paper the same strategy as in Ref. |IJ will be employed to solve the evolution 
equation. Namely, we will solve it in the two kinematical regions where: (i) 7 is small (Q 2 > Q 2 , 
Q 2 ^> A 2 ) and (ii) 7 is close to unity (Q 2 -C Q 2 ). Then the two solutions will be matched in the 
intermediate region. It turns out that the solution inside the saturation region depends strongly 
on the initial parton density inside nucleus. In the Fig. [I] we show two possible situations. If the 
nucleus is light, then the Glauber initial condition is out of the saturation region. We argued 
in Ref. (and will re-derive in this paper again) that the parton density is constant along 
the critical line. This condition sets the initial condition for the solution inside the saturation 
region. As the result, this solution is a function of only one variable. This case was studied in 
Ref. ||1], 0] and will be derived again using another approach in Sec. |3|. On the other hand, if 
the initial parton density is large (because the nucleus is very heavy), then one has to match 
solution inside the saturation region to both the condition on the critical line BC (y > yo), 
and Glauber initial condition AB (y = yo), see Fig. [l|. This cannot be done with the scaling 
solution, so the scaling solution breaks down for very heavy nuclei. We will investigate both 
cases in Sec. [| 

As usually we use the nonlinear evolution equation derived in Ref. in the framework of the 
dipole model |T7|]. Since colour dipoles diagonalize the scattering matrix at high energies |T9|, |T? 



we consider this approach as the most clear from the physical point of view. The evolution 
equation for the dipole-nucleus forward scattering amplitude N{k, y) in the momentum repre- 
sentation readsQ 

F)N(h iA 

= a s x(i(k))N(k,y)-a s N 2 (k,y) , (3) 



dy 

where x (t(&)) is an operator such that 

is an operator corresponding to the anomalous dimension of the gluon structure function and 
the operator x corresponds to the following function 

X (7) = 2^(1) -^(1-7)"^ (7) (5) 



* Eq. ([5|) depends on the impact parameter b t only parametrically via the initial condition Eq. The 
question of bt dependence was discussed in details in Ref. jjj. In this paper we address the case bt = which is 
sufficient for our arguments. 
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which is an eigenvalue of the the BFKL equation. Eq. (^) is derived assuming that the typical 
transverse extent of the dipole amplitude is much smaller then the size of the target and the 
typical impact parameter x <C Ra, h- 

The amplitude in the dipole configuration space is found by employing the Fourier transform 

N(x,y) = x 2 / dkkJ (kx)N(k,y) (6) 
Jo 

N(k,y) = / —J (kx)N(x,y) . (7) 
Jo x 

Initial condition for this equation in the dipole configuration space reads (see Eq. (|])) 

N(x,y = y ) = l-e~ c ^- r+r ^ , (8) 

where r A = In ( 5,Afl y^ )S(0) ) and r = ln(Q 2 /4A 2 ) for Q 2 > Q 2 . In the double logarithmic 
approximation one can neglect the logarithmic dependence of xbG(xb, 4/x 2 ) as well as a s on 
the dipole size, and in addition Q 2 « 4/x 2 . 

Eq. @ can be Fourier transformed to the momentum space where it reads 

1 / k 2 p~ rA \ 

where T is incomplete gamma function of zeroth order. 



2 Solution at Q 2 > Q 2 s {xb)> 

Function x{l) ls singular at its end-points 7 = 0,1. The first one corresponds to the kinematical 
region Q 2 > Q 2 {y), whereas the second to the kinematical region Q 2 Q 2 (y)- Expansion 
of xil) near ^ s end-points and taking the leading contribution corresponds to the double 
logarithmic approximation, i.e. to the summation of a"ln n (<5 2 /A) ln n (l/xs) terms to the right 
of the critical line, and to the summation of a" ln n (Q 2 /Q 2 ) ln n (l/xs) terms to the left of it. 
In the right vicinity of the 7 = one obtains using Eq. (Q) 

11 

X=- = t— r • 10 

7 1 + O r 

Eq. (^) in this approximation reads 

N y (r,y) = a s (l + d r )- 1 N(r,y)-N(r,y) 2 . (11) 

In this paper we will always denote the partial derivatives by subscripts. Note, that the following 
equation 

(1 + d r ) ■ ^ f dk' 2 N(k', y) = N(k, y) (12) 
k z Ja 2 

identifies the operator inverse to 1 + d r . Introducing the new function G by 

G(k,y)= [ k dk' 2 N(k',y) , (13) 

J A 2 
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one gets for it the following equation: 

1 



G ry (r, y) = a s G(r, y) ~ e- r (G r (r, y)) 2 — . (14) 

The advantage of introducing the function G is that we got the nonlinear evolution equation 
in which the nonlinear term has the simplest form. 

Having in mind the success of the semiclassical approach to solution of the DGLAP evolution 
equation at low xb, we look for the solution to the Eq. flli]) in the form 



G(r, y) = e s{r ' y) A 2 , (15) 

where the function S is supposed to be such that \S r S y \ ^> \S ry \. Equation for S(r,y) is 

S r (r,y)S y (r,y) = a a - e' r+s ^S 2 (r,y) . (16) 

This is a partial differential equation of the first order and can be solved by characteristics 
method. Denote S r = f and S y = u>. Then, the characteristics of the Eq. (|16|) can be found by 



solving the system of ordinary differential equations [20 



dr 

~dt 

dy 

dt 
dS 

~dl 
d[_ 

dt 

doj 

~dt 



u + 2fe- r+s , (17) 

/ , (18) 

2fu + 2f 2 e~ r+s , (19) 

e - r+Sf2(1 _ f) ^ m 

-ue- r+s f 2 , (21) 

where t is the parameter along the characteristic. Using Eq. fllED in Eq. ( |19| ) one obtains 

S = 2a s t + S . (22) 

By subscript we distinguish the integration constants. It follows from Eq. (|20| ) and Eq. (|2T| ) 
that 

/ - 1 = uo = Coo . (23) 
u 



Substituting Eq. (^) and Eq. (|23|) into Eq. fl2~T|) results in equation 

— = (Clu 2 + lu — a s )u = Clu(uj — uj + )(uj — ujJ) , (24) 

(aiL 

where 

-1± Vl + ACa s 

«* = 2C • (25) 

Consider some characteristic at large values of r (at t — > oo). In that region nonlinear term 
is small and Eq. (^l|) implies u = u = const. Moreover, f is small and thus u = a s / f is 
large. It is straightforward to calculate u + and u;_ neglecting nonlinear terms 

u> + = ujo and u_ = — s . (26) 

uj - a s 



This shows that in the kinematical region of large r, where nonlinear term is small, solution 
of the system Eqs. ( |T7| - |2"TD can be obtained by setting uj = u = uj + which yields us with the 
usual double logarithmic approximation to the DGLAP and BFKL equations. 

As t decreases uj will no longer be equal to ujq. Corrections stem from the nonlinear term 
(but oj = uj + still holds as it gives correct asymptotic behaviour of solution). Let us address the 
question of what happens if the deviation from the linear behaviour is small (this is equivalent 
to checking the stability of the linear solution). In other words, assume that uj = u + 8u where 
5uj <C uj . Inserting this into the Eq. (p4|) and integrating gives 



u(t) =uj (l + f3e-^- 2 ^ 1 ) , (27) 
where (3 is the integration constant, and by Eq. (|23|) 



f{t) = ^l(l- ^_^l^ e -^-2a 3 )t ) (2g) 

uj \ a s 



Finally, using Eq. (|T|),Eq. (g7|) and Eq. in Eq. ([T^) and Eq. ((TSJ) , the solution to the 
system Eqs. ([T7|-^]]) at large t reads 

r-r = u> t- . 2 "°: 3 ^ /3(e-^-^-l) , (29) 
(uj ~ 2a s )a s v ' 

2/ -2/o = f^ Rie-^-^-l) . (30) 

^o(^o - 2« s ) v ' 

It is readily seen that if t^o > 2a s the solution for characteristic is stable, which means 
that all nonlinear effects which are introduced through the (Glauber) initial condition die out 
exponentially as t increases. By contrast, at ojq < 2a s the solution is not stable. 

This discussion reveals a remarkable property of nonlinear evolution equation: there exist 
such characteristic, that divides the whole kinematical region into two patches[0]. One in which 
all nonlinear effects are negligible, another in which saturation occurs. This line is given by 
condition uo = 2d s . Using Eq. (|29|) and Eq. (|30|) we get the familiar equation for the critical 
line!- 

r - r = Aa s (y - y ) . (31) 

In the next sections we will denote this line as r s (y) or y s {r). 

Excluding u and t (at r — > oo) from Eq. ([2"9"D, Eq. ([ft]) and Eq. Q2"2"D we arrive at the general 
solution of Eq. (|16D in this limit: 

S(r,y) = 2^a s (r-r )(y-y ) + S . (32) 

On the critical line S\ s — r — ro + Sq (we denote by symbol \ s all quantities calculated on the 
critical line). Therefore, by Eq. (|l"3|) and Eq. (|l5|) 



N(k,y)\ s = L Gr {k,y)\ s = e s °- ro = const. (33) 

The dipole scattering amplitude is constant on the critical line. The value of this constant is 
given by Eq. @: 

N(y = y ,r = r A ) = ~T (l)n0.1 , (34) 
in accord with our expectations that the critical line is the boundary of the linear region. 
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3 Solution at Q 2 < Q 2 s (xb)- 



In this section we solve the nonlinear evolution equation inside the saturation region. To this 
end expand x hi the left vicinity of the 7 = 1 using Eq. (||) and Eq. (f|) 

x =T^=~i ' (35) 

This results in the following equation (see Eq. (|3])) 

- N yr (r, y) = a s N(r, y) (l + 2iV r (r, y)) . (36) 
Following Ref. [TJ we introduce the auxiliary function ip by 

N{r,y) = \[ dr' (l-e-^)) . (37) 

Note, that r = \ia(Q 2 /Q 2 ) < is defined with respect to the Q 2 which is the largest scale in 
this kinematical region. For the function ip Eq. (|36|) takes especially simple form 

<Pyr = -a. (l - e~^) . (38) 

The idea of the solution is to split this equation into two and then to find such solution 
of one of the equations that satisfy the other as well. We are looking for the solution in the 
form ip = (p(e s ), where S is assumed to be semiclassical, i.e. |<S ra | |<5 r «S 2/ |. In other words 
the function tp is a quickly varying function of a slowly varying one S. The form of splitting is 
inspired by the solution of the previous section. Namely, we restrict solutions of Eq. (^) by 
the requirement that 

S r (r,y)S y (r,y) = —C = const. (39) 

In other words we assume that e s satisfies the linearized equation (i.e. <p — > 0). 
In the semiclassical approximation defined above, Eq. (^) can be rewritten as 

^ = |(i-^) ■ <*» 

To solve this equation we interchange the independent S and dependent ip variables 0. Denote 
Then, Eq. (|40|) reads 

^CM=|(l-e^>) , (42) 
and can be easily integrated for ( and then for ip. the final result is 

S = I ; y + const. (43) 
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Now we turn to Eq. (|39|) . It is merely the DGLAP equation in the double logarithmic 
semiclassical approximation. Applying the method of characteristics one arrives at the following 
system of equations 

r(t,r) = u' {r)t + r' {r) , y(t,r) = &(r)t + y (r) , S{t,r) = -2tC + S {r) , 

(44) 

f'(t,r) = /o(t) = const. , u'(t,r) = uj' (t) = const. 

where we introduced the prime notation to distinguish variables of this sections from those in 
the previous one. r parameterizes the initial condition which we specify below. By analogy 
with Eq. fl32|), Eq. ( [44] ) give the general solution 



S(r,y) = 2VcJ-(r-r' )(y-y )+S . (45) 



Let us now specify the initial conditions for Eq. (|38|). They depend on the initial parton 
density inside a nucleus. If the initial density is small (curve 1 in the Fig. |l|) then the results 
of Sec. |2] imply that 

<f(y,r)\ s = const. (46) 

The initial condition is different in the case of large initial parton density (curve 2). In this 
region the initial condition reads: 

<p(v,r)\ B = const. , if y Q < y = y s {r) , (47) 

<p(y,r) = -l n (l- e — p(™)) , if y = y . (48) 

Eq. ([4^) follows from its definition Eq. ( p7|) and Eq. (|9]). 

Generally, to satisfy initial conditions one has to find such a surface among those given 
by the solution Eq. (El) which paths through the line specified by Eq. fl4"6] ) or Eq. (B7p and 
Eq. (f48D . To do this one should parameterize the initial condition using parameter r. Then 
calculate /o(r) and oj' {t) using the following equations! 



~d^- fo ^ + UJ °^ ' (49) 

and 

fWo = -C ■ (50) 

The desired solution of the Cauchy problem S(y, r) is then obtained by excluding t and r from 
Eq. (@). 

Eq. ([|^) and Eq. (|5C]) can be easily solved numerically, but to give a transparent analytic 
treatment we notice that the nonlinear corrections in the kinematical region Q 2 < Q 2 S are so 
large that e - ^ 1. Indeed, by Eq. (^) and Eq. (0), the value of ip on the critical line 

r = r s (y) is 

ip Q = -ln(l-e _1 ) ^0.5 , (51) 

which follows from the requirement that the point (ta, yo) lies on the saturation scale (see 
Fig. 0). Then Eq. ( p8|) can be easily integrated 

¥(y,r) = -a s yr + f^y) + f 2 (r) , (52) 

where fi and f 2 are arbitrary functions. 
Consider two different initial conditions. 



S 



1. In the case of large initial parton density the initial condition Eqs. 

<p(r,y = y ) = -a s y r + h{ m ) + f 2 (r) = - In (l - e - exp ^) , (53) 
<p(r,y = y,(r)) = -a s y s (r) r + fi(y s (r)) + / 2 (r) = <p , (54) 

and is satisfied by the following function 

(p{y,r) = -a s {y-y )r + a s {y-y )r s {y)-ln(l-e- exp ^ rA ^ + 

\n(l-e- exp ^ y) - rA) )+cp . (55) 

2. If the initial parton density is small, then in addition to the initial condition Eq. ( |4q) for 
the function (p, we have to specify the initial condition for its derivative. Note, that to 
the right of the critical line p = 2^Ja s (y — y )r — r + (p . Thus ip r \ s = — |. This yields 

<p(r,y = y s (r)) = -a s y s (r) r + fi(y s (r)) + f 2 (r) = <p , (56) 

~2 = -®sy s (r) + /ar(r) . (57) 

The solution is 

<p{r,y) = i(r - ^a s y f - -(r - ^a s y) + <p = ^z 2 - + (p . (58) 

It coincides with the asymptotic solution obtained in Ref. 

We see how the introduction of the auxiliary function (p simplifies the solution of Eq. (|36|) . 
The main point is that while the value of the amplitude iV on the critical lineEq. (j34|) is 
intermediate between perturbative and saturation regions and thus all terms in Eq. (|36"D are of 
the same order, the value of (p on the same line Eq. (|5T|) is such that Eq. (^) can be greatly 
simplified: all nonlinear corrections to iV are described with good approximation by just a 
constant in the equation for p. 



47, 4|) read 



4 Discussion 



In Sec. 2 we found that in the kinematical region to the right of the critical line the nonlinear 
corrections are small in the double logarithmic approximation. This result is in agreement with 
Ref. [|H]. The evolution equation was solved in the semiclassical approach using the method of 
characteristics. The characteristics intersect the line y = yo with different angles, the tangents 
of which are given by the values of fo- At low /o (large c^o) characteristics run in the linear 
kinematical region. As ujq decreases characteristics approach the critical line and finally, at 
uj = 2d s there exists the characteristic which runs along the boundary of the saturation region. 
In the double logarithmic approximation the characteristic with u = 2a s is merely the critical 
line. On this critical line the scattering amplitude was found to be constant (see Eq. (|33|)). 

To find the A-dependence of the saturation scale we have to require that the critical line 
(which is the characteristic at large r) satisfy the following condition 

r s (y = Vo) = r A , (59) 
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which means that the Glauber initial condition Eq. ([8]) itself provides a saturation scale (r^). 
However, when r s ~ the value of uq may differ from that at r — ► oo (i.e. 2a s ). This is 
because our approximation that \S r S y \ ^> \S ry \ which led to Eq. (|T6| ) breaks down at r s « r^. 



The exact solution of the (linearized) evolution equation which satisfies the (linearized) initial 
condition (see Eq. (|8|)) is 

jf^yj = J ffi c -T-+r A +a s (i/-i/o)^+T-/ DLA c 2^/g s rfa- ? ;o)+r A -r ^qq^ 

(the integartion contour runs along the line parallel to the imaginary axes to the right of all 
singularities of the integrand). Assuming that the remarkable property N\ s = const, holds up 
to r = ta one gets (neglecting logarithmic contributions) 



2^a s r(y - y ) + r A - r = <p = const. (61) 

Thus, 

a s (y s (r) - yo) = {r 'l A? , (62) 

where = to satisfy r s (y Q ) = r&. At r ^> the A dependance of the saturation scale is 
given by 

r = Aa s y + 2r A , (63) 

in full correspondance with our previous workfllj and esimations in the Introduction. 

In Sec. 3 we showed that the initial parton density in the nucleus plays very important role 
in the saturation regeme dynamics. Namely, if A is such that initial parton density is smaller 
then some critical value, then the solution inside the saturation region depends only on one 
variable Q 2 /Q 2 s (xb)- This scaling behaviour was predicted by many authors 0, §, [14], [Tj| and 



observed in experimental data|[L8|l. It turns out that if the initial parton density is larger then 
the critical value, the scaling behaviour breaks down. Let us estimate at what A the parton 
density reaches its critical value. It is seen in Fig. [I] that this occurs when = r , where r 
defines the initial point of evolution, i.e. it satisfies a s (r ) <^ 1. Using definition of we get 

^dipoleS(O) = 1 , (64) 

where Odipoie is a cross section for the small-dipole-nucleon interaction. Since at xb = 10~ 2 (i.e. 
V = Vo) Odipoie ~ Icvrpl^T]] and S(0) = 2RaP, where p is a nuclear density, we have 



^a 7Tp pR A = l , (65) 

Taking a wp = 24 mb, p = 0.17 fm~ 3 and Ra = l.lA 1 ^ 3 fm one gets A cr it. ~ 70. The direct 
estimation using double logarithmic formulae are not reliable since they are valid up to some 
logarithmic corrections and perhaps a numerical factor, which are quite important since A 
enters to the 1/3 power. 

Concluding this paper, we would like to emphasize that we found that the saturation scale 
for DIS on nuclei lighter then A ~ 70 scales like Q 2 S ~ A 2 / 3 at all energies, whereas for DIS on 
nuclei heavier then A ~ 70 it scales like Q 2 ~ A 1 / 3 when the evolution starts and then reaches 
the same scaling behaviour as for light nuclei Q 2 ~ A 2 / 3 at very high energy. 
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